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The theory of quantum error correction provides a general methodology for protecting quantum 
^"^ information from noise. It is therefore expected that quantum error correction will be essential 

^ in operating future quantum computers, machines storing and manipulating very large amounts of 

CN quantum information in the course of long quantum computations. 

Certainly, there is nothing particularly quantum mechanical in the idea of encoding the infor- 
mation stored and processed inside computers. Ordinary digital computers already use various/aw/^ 
VO tolerance methods at the software level to correct errors during the storage or the transmission of 

information — e.g., the integrity of the bits stored in hard disks is verified by using parity checks 
r—i (checksums). In addition, for critical computing systems such as these inside airplanes or nuclear 

O reactors, software fault tolerance methods are also applied during the processing of information — 

e.g., airplane control computers compare the results from multiple parallel processors to detect 
faults. In general however, the hardware of modern digital computers is remarkably robust to noise 
^ so that, for most applications, the use of additional software error correction is rather limited. 

, ^ In contrast to the easiness and robustness with which classical information can be processecQ 

the processing of quantum information appears at present to be much more challenging. Although 
^ constructing reliable quantum computing hardware is certainly a daunting task, we have neverthe- 

00 less strong hopes that large-scale quantum computers, able to implement useful long computations, 

"^ can in fact be realized. This optimism is founded on methods of quantum fault tolerance which 

^ show that scalable quantum computation is, in principle, possible against a variety of noise pro- 

1*-^ cesses. Demonstrating that these methods work effectively in practice is a major challenge for con- 

O temporary science, a challenge whose outcome will depend on our progress in understanding the 

I physical noise processes in experiments, and on our ability to design and optimize fault tolerance 

1^ methods according to the limitations and the noise characteristics of experimental devices. 

This chapter is an introduction to software methods of quantum fault tolerance. Broadly speak- 
ing, these methods describe strategies for using the noisy hardware components of a quantum com- 
^ puter to perform computations while continually monitoring and actively correcting the hardware 

faults. The methods we will discuss are general and apply independently of how the hardware com- 
ponents are physically realized in the laboratory. Nevertheless, one should not lose sight of the fact 
that what we describe in this chapter as elementary hardware components are not elementary from 
an experimental point of view. Already at the level of the realization of qubits in the laboratory, 
the experimenter strives to choose implementations with high inherent robustness to noise such as 
qubits encoded in decoherence-free subspaces or noiseless subsystems, or qubits which are topo- 
logically protected. In addition, noise in the elementary hardware operations can be suppressed 
by using various open-loop techniques such as refocusing or dynamical decoupling. Even though 

^ Of course, this was not always the case; photographs of ENIAC, the first universal electrical computer, speak volumes about how 
difficult the first steps of classical computing were. 
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these various qubit encodings and noise-suppression techniques can be highly effective, some resid- 
ual noise will always remain; it is this residual effective noise that needs to be treated by the error 
correction and fault tolerance methods we will discuss in this chapter. 

The basic conceptual ideas of fault tolerance for quantum computation are very similar as in 
the case of classical computation: First, a code is chosen and each logical step of the computation 
is implemented by a fault-tolerant gadget which acts on the encoded information; these gadgets 
comprise many elementary hardware operations, and they are designed to implement the desired 
logical transformation on the encoded information while at the same time detecting and correcting 
errors. And secondly, the protection from noise is increased by designing a hierarchy of encoding 
layers such that errors become progressively weaker as we pass from one layer to the next. 

Despite these similarities, there are two major differences between quantum and classical fault 
tolerance, which are related to the differences between classical and quantum error correction. The 
first difference is that in the quantum case error correction needs to be implemented coherently, i.e., 
in a way that preserves the quantum superpositions in the encoded information that is processed by 
the quantum computer — this requirement has no analogue in classical fault tolerance since quantum 
superpositions and quantum interference play no role in classical computation. The second differ- 
ence relates to the types of noise that are of concern in the two cases. For ordinary computers that 
manipulate classical information digitized in bits, noise can simply be viewed as causing abrupt 
changes in the value of each bit (bit flips). For quantum computers on the other hand, information is 

stored in quantum states which (if pure) are in general superpositions aie^^^ |^i) + Ci2e^^'^ |^2) H 

of various physically relevant basis states 1^^) with real coefficients a^ and 0^; then not only can 
noise cause changes in the amplitudes a^ (which is analogous to the bit flip errors for classical in- 
formation as af is the probability of occupation of the state 1^^)) but noise can also cause changes 
to the phases (pi (phases are irrelevant when storing classical information but they are important 
quantum mechanically as they determine the ability of the superposed basis states to interfere). 

1.1 Quantum circuits and error discretization 

The precise character of noise in future quantum computers will depend on the particular hardware 
implementation. There is a wide variety of prospective implementation schemes that are being ex- 
perimentally investigated at present, but in this chapter we restrict the discussion to those schemes 
which fall under the quantum circuit model for which methods of fault tolerance are better under- 
stooqB 

Quantum circuits are a generalization of classical circuits: A classical circuit computing a 
boolean function f on n bits is a prescription for expressing / as a composition of functions or 
gates that act on a fixed, independent of n, number of bits at a time; gates are chosen from a finite 
set which is universal allowing any function to be computed — e.g., the NOT gate which flips the 
value of a bit together with the AND gate which computes the conjunction of the value of two bits 
form a universal gate set. Similarly, a quantum circuit is a prescription for implementing a physical 
operation on the Hilbert space T-Ln = (C^)*^^ of n qubits as a composition of elementary physical 
operations which are applied on a fixed number of qubits at a time. Although Hn is continuous, there 
exist finite sets of physical operations acting on at most two qubits which are quantum universal al- 
lowing the approximation of any physical operation in T-Ln to any desired accuracy; these universal 
sets comprise preparations of single qubits in certain pure states, certain unitary transformations 
or quantum gates on single qubits or between pairs of qubits, and measurements of single-qubit 
observables. Our diagranmiatic representation of a quantum circuit is shown in fig. |l.l[ 

^ In particular, we will not discuss quantum computations realized purely by adiabatic evolution for which a general theory of fault 
tolerance is lacking. 
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Fig. 1.1. A quantum circuit is a sequence of qubit preparations (7^), qubit measurements (A^), and quantum 
gates (a • denotes a gate applied on a single qubit, and two •'s connected vertically denote a gate applied on 
a pair of qubits). The measurement outcomes are processed by classical computers alongside the quantum 
computer; intermediate measurement outcomes condition the application of future quantum gates, while the 
final measurement outcomes encode the answer of the quantum computation (0 or 1). 



We will not discuss examples of quantum universal sets here. Conceptually, what is important is 
that a finite number of elementary operations suffices for implementing any quantum computation. 
Therefore, a quantum computer is a discrete machine just like a classical digital computer; in the 
classical case, the elementary hardware components are gates on one or a few bits, while in the 
quantum case they are physical operations on one or a few qubits. Discreteness is essential both 
classically and quantumly because it implies that fault tolerance can be achieved by constructing 
fault-tolerant gadgets for each one of the operations in a universal set and composing these gadgets 
together. 

The second essential ingredient for fault tolerance is the ability to discretize errors so that error 
correction becomes possible. In classical digital computation, the basis of error discretization is 
the digital encoding of information. Once a bit of information is represented in a physical quantity 
taking the value v^ to encode and v\ to encode 1, noise in gates can be described in terms of 
discrete errors taking v^ to v\ or vice versa, and small fluctuations around the values v^ and v\ 
can in practice be ignored — e.g., a NOT gate can be implemented by a CMOS inverter in saturation; 
the input and output bit values and 1 are encoded as different voltages v^ and v\, and the output 
voltage is essentially insensitive to small variations 8v <^v^^ v\ in the input voltage. 

In the case of quantum computation, it is important to recognize that there is no unambigu- 
ous way to say which qubits of an entangled multi-qubit quantum state processed by the quantum 
computer are erroneous and which are not — e.g., consider the maximally entangled two-qubit state 

|$o) = ^(|0)|0) + |l)|l)); (1.1) 

if noise acts on |$o). there is no good way to say which of the two qubits is erroneous because for 
any single-qubit operator E, {E ® /)|$o) = (/ ® £^^)|$o), and whether the state has suffered an 
error cannot be determined by just observing the properties of any one of the two qubits in isolation. 
It is therefore helpful to avoid using a semantic language where the notion of an error depends on 
the quantum state on which errors act; instead, we adopt a syntactic language where the notion of an 
error is defined operationally independent of the actual quantum state. Instead of associating errors 
with individual qubits, we can associate errors with mutually orthogonal subspaces of the entire 
Hilbert space of a collection or a block of several encoded qubits; for a block comprising n encoded 
qubits 

T~Ln — ^J7^n 5 (1-2) 

s 

where the superscript 5, called the syndrome, is a label for the different subspaces 7/^ for the n- 



qubit Hilbert space Hw One of the subspaces, the code space T-L^, is the preferred one in the 
sense that quantum information is encoded in a quantum state that is supported in H^. Because of 
noise, the encoded quantum state will tend to escape from the code space toward other subspaces. 
But because all subspaces are mutually orthogonal, there is a generalized measurement, called a 
syndrome measurement, which allows different subspaces to be distinguished unambiguously. By 
performing a syndrome measurement, we can then use the measurement outcome fi to deduce s and 
therefore learn about the subspace on which the noisy encoded quantum state is supported. Since 
there is only a finite number of orthogonal subspaces, we can execute the computation that yields 
s given ji digitally by processing the measurement outcome in a classical computer. If the result 
of this computation is a nontrivial syndrome value which indicates that the noisy encoded quantum 
state is supported in a subspace different than the code space, we can apply a recovery operation on 
the encoded qubits, which is conditioned on the result of the classical computation and which returns 
the encoded quantum state to the code space. Because classical digital computers are in practice 
extremely robust to noise, we usually assume that the classical processing of fi to obtain s and the 
classical control of the quantum computer conditioned on s can be implemented perfectly without 
faults. Of course, no matter how improbable, faults in the classical processing of the measurement 
outcomes can lead to a failure to apply the appropriate recovery operation, so that the accuracy of 
the quantum computer is ultimately limited by the accuracy of the on-the-side classical computer. 



1.2 Noisy quantum computers 

We have seen that discretizing the entire Hilbert space of the qubits processed by the quantum 
computer allows us to encode quantum information in a quantum state supported in the code space 
and to protect against noise that takes this encoded quantum state to other orthogonal subspaces. 
But for what types of noise processes is this method of encoding quantum information effective? 
Certainly, we have very little hope of protecting information against noise that acts collectively 
on many hardware components of the quantum computer and whose strength is not moderated as 
a function of the number of qubits it affects — e.g., we would be helpless if a power outage, an 
earthquake, or a high-energy cosmic ra}Qwere to hit our quantum computer affecting many qubits 
all at once. 

If we exclude such malicious types of collective noise against which no error-correction method 
can be effective neither for quantum nor for classical computation, we are left with several other 
contributions to noise that need to be considered: First, there is noise due to imprecisions in the 
implementation of each elementary hardware operation — e.g., noise in the control parameters dur- 
ing the implementation of a unitary gate U might result in realizing another unitary U+5U instead, 
where 5U may be systematic or it may vary stochastically. Secondly, there are unwanted inter- 
actions among the qubits in the quantum computer — e.g., an electromagnetic coupling of nearby 
quantum-dot qubits which decays as a power of their relative distance. Thirdly, there are interac- 
tions between the qubits of the quantum computer and an environment represeniing external degrees 
of freedom which are not under our control — e.g., a coupling of integrated superconducting qubits 
to nuclear spins in the substrate. And finally, in settings where qubits are realized by selecting 
a two-dimensional subspace inside a multi-dimensional system, there is noise which couples the 
two-level qubit subspace to other levels of the same system — e.g., if a qubit is realized by using 
the ground-state hyperfine splittings of a trapped ion, noise can induce transitions between these 
hyperfine levels and other higher-energy levels of the ion. 

A useful classification of these noise processes concerns their spatial and temporal locality. 

^ High-energy cosmic rates are unlikely on earth, but they are a real concern for computers inside space shuttles. 



1.2 Noisy quantum computers 



Intuitively, we say that noise is spatially local or simply local if, during any time interval, it acts col- 
lectively only on qubits which are interacting in the ideal quantum circuit at the same time interval — 
i.e., if two-qubit gates are applied in parallel to several pairs of qubits during a specific time interval, 
local noise can act collectively on qubits that belong to the same pair but not on those that belong 
to different pairs. Locality is a desirable property because it implies that noise cannot afflict global 
damage by causing the simultaneous failure of many hardware components. Although in general 
error correction fails for non-local noise, there are in fact some non-local noise processes for which 
effective error correction is possible: First, certain types of collective noise whose nature is known 
in advance can be suppressed effectively by using the techniques of decoherence-free subspaces and 
noiseless subsystems. And secondly, there are types of non-local noise which can be treated as if 
noise were local and for which fault-tolerance methods designed to protect against local noise are 



effective; we will discuss two such examples in Section [L23 

With regard to temporal locality, the question is how correlated is the noise that acts on different 
hardware components which are executed at different time intervals. We say that noise is temporally 
local or Markovian if the noisy evolution can be described by using a sequence of superoperators 
each taking the density matrix at the end of one time interval to the density matrix at the end of the 
following time interval. 

When the quantum computer implements a quantum circuit of depth[^i?, we can discretize the 
total computation time T inD intervals Ti, T2, . . . , Td each of duration to equal to the time it takes 
to execute an elementary operation. The Markovian property then translates to the requirement that 
noise has a typical correlation time comparable to to. Since the interaction between the quantum 
computer and the environment is incoherent across different time intervals, we can trace over the 
state of the environmental degrees of freedom after each interval to obtain a reduced density matrix 
describing the state of the quantum computer. In this case, the noisy evolution is described as a 
mapping between the reduced density matrices at different intervals; 

Pj ^Sj{pj-i) , (1.3) 

where pj is the reduced density matrix at the end of interval Tj, and Sj is a superoperator describing 
the evolution from the end of interval Tj_i to the end of interval Tj. 

On the other hand, if there are noise processes with typical correlation times longer than to, 
we cannot obtain an accurate description of the noisy evolution by tracing out the external degrees 
of freedom after every time interval. In this case, the information that the environment exchanges 
with the quantum computer could in principle be retained for long times so that we cannot simply 
describe the entire noisy evolution as a composition of superoperators. Because of this reason, the 
analysis of the effects of non-Markovian noise is more demanding than for simple Markovian noise 
and, as we shall see in the next section, our conclusions for the effectiveness of fault-tolerance 
methods against non-Markovian noise are generally weaker than for Markovian noise. 

In this section, we will discuss several concrete examples of noise models that have been ana- 
lyzed in the context of fault-tolerant quantum computation. 



1,2,1 Setup 

We consider the noisy implementation of an ideal quantum circuit comprising L elementary opera- 
tions followed by the final qubit measurements whose outcome encodes the result of the computa- 

^ The depth of a quantum circuit is the maximum number of elementary operations appHed on any qubit of the quantum computer 
(including the identity operation which is implicitly applied when a qubit is stored while operations are applied on other qubits). 



tion. The ideal quantum circuit produces the quantum state 

f}^^^^^OLO---o02oOi, (1.4) 

where the superoperators Oj correspond to either a qubit preparation, a unitary gate, or an interme- 
diate qubit measurement which conditions subsequent operations; the operation for the preparation 
of a qubit in the pure state |^) is 

V\^) = |V)(V'I , (1.5) 

the operation for a quantum gate U appHed on input X is 

U{X) = UXU^ , (1.6) 

and the operation for a projective measurement of an observable a applied on input X with projector 
Ma corresponding to measurement outcome a, 

a 

For simplicity and since noise on the final measurements can be modeled by noise acting in the 
immediately preceding operations, we may assume that the final measurements are implemented 
ideally without faults. Finally, as we have noted, we will consider performing the processing of 
the outcomes of both the intermediate and the final measurements in a classical computer operating 
alongside the quantum computer, and we will assume there are practically no faults in this classical 
hardware. 

7.2.2 Local Markovian noise 

Our first example is noise which is both local and Markovian. The Markovian property implies that 
we can describe the noisy evolution as a composition of superoperators, and the locality property 
implies that the superoperator for each time interval can be expressed as a tensor product of super- 
operators, each superoperator corresponding to one of the different elementary operations which are 
implemented in parallel during that interval. 

We may express the superoperator describing the noisy implementation of each elementary 
operation as J\fj o Oj, where Oj is the ideal superoperator and Mj is a superoperator describing 
deviations from the ideal due to noise — by definition, the support of J\fj is contained in the support 
of Oj when noise is locaF] Therefore, because of noise, instead of the ideal quantum state in 



eq. ( |1.4| ), the quantum computer really prepares the state 

^noisy = AAl o O^ o . . . o AA2 o O2 oM o Oi , (1.8) 

where o denotes composition. 

Since all J\fj would be trivial were there no noise, a natural measure for the noise is the distance 
between J\fj and the identity superoperator X, and we can define the noise strength 

^^maxllAAj -Xllo , (1.9) 

3 

where 1 1 • | |o is a suitable superoperator norm, the diamond nornr] If we now write Mj^X^ Tj 

^ The support of a superoperator O (or operator O) acting on density matrices (or quantum states respectively) defined on various 
subsystems is the tensor product of the Hilbert spaces of all the subsystems on which O (or O respectively) acts nontrivially; in our 
case, the subsystems are the qubits of the quantum computer and any subsystems in the environment. 

■^ If the superoperator E has an n-qubit input, and Xn is the identity superoperator on n qubits, then 1 1 ^ 1 1 o = max 1 1 (^ (g) /n ) (X) 1 1 1 , 
where we maximize over all X such that | |X| | i=TrvXtX=l. Note that we could optimize our estimate for the noise strength by 
taking X to be proportional to the identity superoperator with a proportionality constant of magnitude between and 1 ; but here we 
will not discuss this generalization. 



1.2 Noisy quantum computers 7 

for some fault operator, or simply fault, Tj and substitute in eq. ( |1.8[ ), we find 

^noisy ^ (X+ J-^) o Ol o . . . o (X + 7-2) o O2 o (X+ J-i) o Oi . (1.10) 

By opening all parentheses, we obtain a sum of terms corresponding to diffcreni fault paths; in each 
fault path, faults Tj have occurred in a specific subset of the L elementary operations, while the 
identity superoperators are applied on all remaining operations. In particular, we can write 

^noisy _ ^ideal ^ ^faulty ^ ^j jj^ 

where p^^^^^ corresponds to the unique fault path where identity superoperators are applied every- 
where and (^f^^i^y contains all other fault paths for which there is at least one insertion of a fault — we 
use ( instead of p in (^^^^^^y to emphasize that it is not a density matrix but rather the difference of 
two density matrices. 

The fault path expansion is helpful for understanding how accurate is the noisy circuit. More 
precisely, we would like to know what is the distance 6 between the probability distribution {q^^^^^} 
for the outcomes {/i} of the final measurements on the noisy circuit and the distribution {gjf^^^} if 
these measurements were applied on the ideal circuit instead. We can express S in terms of the 
1-norm (or Kolmogorov distance) between the two probability distributions, 

S = E IC"' - i"*"'! = E |Tr {m, (p--y - p'd-1)) I , (1.12) 

/i /i 

where the projectors M^ define the final measurementaM But next, S can be related to the trace 
norm between p^^^^y and p^^^^^ ; i.e., 

^ = EiE^-<^-i^/^i^-)i ^ EEi^-ii<^-i^-i^-)i ^ Ei^-i = iic'""'*''iii , (1-13) 

where 1^;^^) are the eigenvectors of i^^^^^^y = ^noisy _ ^ideai ^j|-]^ corresponding eigenvalues Vy. 

Now, how do we upper bound the norm of ^^^^i^y ? if we let C denote the set of all L elementary 
operations in the quantum circuit, 

L 
^faulty ^ ^(_1)-1 ^ C{Cr) , (1.14) 

r=l CrCC 

where the second sum is over all subsets Cr of C with cardinality r, and C{Cr) denotes a sum of 
all the fault paths with faults applied on all operations in the set Cr. Eq. ( |1.14| ) can be derived from 
the inclusion-exclusion trick of combinatorics: Since (^^^^^^y is the sum of all the fault paths with 
at least one fault, the sum of all C{Ci) counts correctly all the fault paths with exactly one fault 
but overcounts the fault paths with at least two faults; to amend the overcounting, we subtract the 
sum of all C{C2) which corrects the overcounting of all the fault paths with exactly two faults but 
introduces an undercounting of the fault paths with at least three faults; and so on. 

For each specific set Cr, C{Cr) is nothing by the composition of the ideal superoperators Oj 
interspersed with faults Tj applied on all operations in Cr and the full noise superoperators J\fj 
applied on all the remaining operations. Of course, since superoperators have unity norm, | |Oj | |o = 
||A/^-||o = 1 and thus 

\\aCr)\\i<e' . (1.15) 



By using the triangle inequality, and since there are (^) distinct subsets Cr of C, eqs. (1.13) and 

^ The projectors M^ are non-negative (i.e., {u\M^\u) > for any \u)), and they are normahzed so that Eu ^m — ^- 



( |1.14[ ) now imply that 

L 



where in the last two steps we assumed that e < 1/L. 



Our derivation of eq. ( |1.16| ) via eq. ( |1.14| ) was made in order to introduce a simple application of 
the inclusion-exclusion trick that is also being used later in this chapter. An improved upper bound 
on 6 can in fact be derived without the assumption e < 1/L by simply noting that we can group 
fault paths depending on their earliest faulty operation; 

L 
^i.nlty^J2aOr). (1.17) 

r=l 

where C{Or) is the composition of the ideal superoperators Oj interspersed with the identity super- 
operators applied on the operations 1 to r— 1, a fault J> applied on the r-th operation, and the full 
noise superoperators J\fj applied on the operations r+1 to L. Since superoperators have unity norm, 



||C(Or)||o < ^ and thus eq. (1.17) implies 



6<L£. (1.18) 

We conclude that, for a constant error strength e, the accuracy 1 — S of the noisy quantum circuit 
decreases at most linearly with the circuit size L, in accordance to what is expected for a discrete 
model of computatioij^ Of course, as we shall discuss in the following sections, the goal of imple- 
menting the quantum computation by using fault tolerance methods is to replace e in eq. ( |1.16| ) by 
a smaller — in fact, an arbitrarily small — effective noise strength, thus making the accuracy of the 
noisy circuit approach as close to unity as desired. 

1.2.2.1 Assessment and examples 

At this point, we can step back to note the two essential assumptions that allowed us to derive 
eq. ( |1.18| ): First, we assumed that the superoperators describing the noisy evolution can be expanded 
perturbatively as a sum over fault paths. Secondly, we assumed that fault paths with many faults are 
exponentially suppressed in the sense of eq. ( |1.15| ). In fact, we may view eq. \\.\5) as the defining 



property of local Markovian noise, even if the superoperators that describe the noisy evolution are 
not strictly local. Thus, we generally say that 

Definition (Local Markovian noise). Noise is local and Markovian if the noisy evolution can 
be expanded as a sum over fault paths, where faults are described as (differences of) superoperators 
and the norm of the sum of all the fault paths with faults in any r specific elementary operations is 
upper bounded by e^ for some constant noise strength e. 

This relaxed definition has the advantage that it can describe correlated noise both in space 
and in time: Subject to the constraint that fault paths must satisfy eq. ( |1.15| ), the fault operators 



comprising each fault path are otherwise unconstrained; in particular, the various fault operators 
can be controlled by an adversary who may chose to act collectively on all the faulty operations any 
way she pleases. 

^ For discrete models of computation, to achieve a constant accuracy 1 — 5, the number of bits of precision required to specify the 
physical parameters associated with each elementary operation to within £ — e.g., the amplitude and timing of a voltage pulse used to 
control a CMOS gate — grows logarithmically with the size L of the computation. (With the bits of precision growing logarithmically 
with L, e decreases polynomially with L.) In contrast, for analog models, the number of bits of precision grows polynomially or 
even exponentially with L. 



1.2 Noisy quantum computers 



The local Markovian noise model captures several noise processes of interest; below, we discuss 
three simple but important examples. 

Control noise. A common source of noise is due to imprecision in the control parameters during the 
implementation of each elementary operation — e.g., noise in the timing or the intensity of external 
magnetic fields used to manipulate the state of a superconducting qubit. 

In the simplest case, consider the implementation of a single-qubit gate corresponding to a 
rotation by an angle 26 around the z direction; this operation is described by the superoperator 

ni{X) = e^^^^Xe-^^^^ . (1.19) 

Because of imprecisions in the control parameters, a rotation by a different angle 29' = 2{9+59) 
for some small fixed deviation 59 may be implemented instead. We can express the noisy superop- 
erator as TZ^, = AActri o T^l^ where AActri = 2^ + J'ctri and 

:F,tr\{X) = iS9(a,X-X a,) + 0{59^) , (1.20) 

so that control noise satisfies eq. \\.\5) with e = 0{59). 

A similar conclusion also holds if the deviation angle is not fixed but varies stochastically, and 
also for control errors in multi-qubit gates, preparations or measurements. 

Relaxation. Another common source of noise is due to thermal relaxation — e.g., in systems where 
|0) and |1) are encoded in different energy eigenlevels, the state |1) may spontaneously relax to the 
lower-energy state |0). To first approximation, relaxation can be expected to act independently on 
each qubit during the execution of a quantum computation; then for each qubit and during each time 
interval, relaxation with a characteristic time scale Ti can be modeled by the amplitude damping 
superoperator AAreiaxl^) = M^XmI + MiXM\, where 

Mo = i±^I+i^4^.. ,M, = ^a.(l-a.), (1.21) 

and 7 = 1 — e~^^l^^ . We can write AAreiax = 2^ + ^reiax where 

^reiax(X) = ^a.XM^ + ^-M^Xc, + M^XM\ + 0(7') , (1.22) 



so that relaxation noise satisfies eq. ( 1.15 ) with 5 = 0(7). 



Probabilistic noise. In many cases noise can be modeled as a random processes — e.g., shot noise 
in the laser fields used to control trapped ionic qubits. Ignoring the details of the underlying random 
process, and letting p denote the probability of a fault during the implementation of each elementary 
operation — if r is the fault rate, p = rto — , noise can be modeled by the superoperator 

AArand(X) = {X - p)X ^ pEXE^^ , (1.23) 

where E is an arbitrary operator acting on the support of the ideal operation (subject to the constraint 
E'^E = / required for AArand to be trace preserving). 
We can write AArand = 2^ + ^rand where 

^rand(X) = -pX + pEXE^^ , (1.24) 

SO that probabilistic noise satisfies eq. ( |1.15 ) with 5 = 0(p). 
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Fig. 1.2. On the left, a unitary gate /7 in a quantum circuit is conditioned on the outcome of a preceding 
single-qubit measurement along the orthonormal basis {|0), |i)}. On the right, a mathematically equivalent 
circuit where the two qubits interact via a unitary gate U followed by a single-qubit measurement along the 
basis {|0), |1)}; here, U = A{U){V ^ I) where the single-qubit unitary V acts on the measured qubit as 
V\b) = |0), V\i) = |1), and the two-qubit A{U) applies U on the second qubit conditioned on the state of 
the measured qubit being 1 1) . 



1.2.3 Local non-Markovian noise 

Our local Markovian noise model is powerful enough to capture several important noise processes 
such as systematic control errors or thermal relaxation. However, requiring that fault paths are 
associated with superoperators is rather limiting because it constrains the possible noise correlations 
between different fault paths; while the fault operators in any specific fault path may be arbitrarily 
correlated, there can only be classical but no quantum correlations — i.e., no quantum interference — 
between the fault operators in different fault paths. 

To go beyond Markovian noise, we can no longer trace over the environmental degrees of free- 
dom to obtain a superoperator description of the noisy evolution. Now, our description will need to 
include explicitly the quantum state of the environment and its joint unitary evolution with the qubits 
of the quantum computer during the course of the entire quantum computation — of course, while 
we assume that we have control over all the qubits of our quantum computer, the environmental 
degrees of freedom are inaccessible and in many cases their precise nature is unknown. 

We assume that the qubits of the quantum computer can be initialized in a pure state |^o)qc — 
e.g., we can prepare all qubits in their lowest-energy eigenstate (at least, come very close to it) by 
cooling — so that the state at the beginning of the quantum computation, including the environment, 
is |^o)qc ® |0o)^ for some unspecified pure stateF]|(/)o)^. If there were no noise, implementing the 
ideal quantum circuit would then correspond to implementing a sequence of unitary operators Uj 
producing the final state 

Wqc' = {Ul---U2- Ui) |Vo)qc , (1-25) 

on which state we finally apply measurements that give the result of the computation. (If there are 
unitary gates that are conditioned on the outcome of a preceding measurement — cf . fig. |1.1] — , we 
can mathematically replace them in our analysis by different unitary gates which are followed by 



measurements as in fig. 1.2 ) 



We may describe the noisy implementation of each elementary unitary operator Uj as Nj • Uj 
where Nj is a unitary operator on the quantum computer and the environment and it describes 
deviations from the ideal due to noise — by definition, the support of Nj is contained in the union of 
the support of Uj and the environment when noise is local. Since all Nj would be trivial were there 

^ We can always obtain a representation of the environment in terms of a pure state since any mixed state can be purified by introducing 
an auxiliary Hilbert space. 
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no noise, we can now define the noise strength as 

e = maxmin||7Vj - Iqc\\oo , (1.26) 

3 Iqc 

where we vary over all unitary Iqc that act trivially on the quantum computer and in an arbitrary 
way on the environment. 



By expanding each Nj as the sum of the Iqc that minimizes the norm in eq. ( 1.26 ) and a fault. 



Nj = Iqc + Fj, we can substitute in eq. ( 1.25 ) to find that the joint state of the quantum computer 



and the environment prior to the final measurements is 

IV')qcI = (Iqc + Fl)-Ul--- {Iqc + F2) ■ U2 ■ {Iqc + Fi) ■ f/i(|^o)QC ® l^o)^) , (1-27) 

where it is understood that the Uj did trivially on the environment. We may now open all the 
parentheses to obtain a fault-path expansion just like we did in the previous section; here, each fault 
path identifies a specific subset of the L elementary operations where faults Fj have occurred, where 
the Fj are differences of unitary operators instead of differences of superoperators as in the case of 
Markovian noise. We can then write 

where for succinctness we have dropped the QCE subscripts but it is understood that all states are 
supported in the Hilbert space of the qubits of the quantum computer and also the environment. 
Here, \^y'^^^^ corresponds to the unique fault path where operators Iqc are applied everywhere and 
I ^^ faulty contains all other fault paths for which there is at least one insertion of a fault — we use i9 
instead of ^ in \'dY^^^^y to emphasize that it is not normalized but it is rather the difference of two 
normalized pure states. 

How can we estimate the accuracy of the noisy quantum circuit in the presence of local non- 
Markovian noise? From eq. ( |1.16| ), it suffices to evaluate the trace norm of the difference between 



the noisy and the ideal density matrices; since the final states are pure, 

5 = iii^)(^r^^'^ - i^)(^r^'^^iii < 2|p)^^^i*y|| , (1.29) 



where we used that |||n)(n| — |z)(z|||i = 2^1 — |(n|z)p < 2|||n) — |z)|| for any normalized pure 
states |n) and \i). It remains to obtain an upper bound on the norm of \i^Y^^^^y , 



In analogy to eq. (1.14), 



L 
l^^faulty ^ J^(_l)-1 Y^ |^(C,)) , (1.30) 

r=l Cr^C 

where \d{Cr)) denotes a sum of all the fault paths with faults applied on all r operations in the 
set Cr\ since |i9(Cr)) is obtained by applying unitary operators everywhere except at the faulty 
operations, 

W{Cr))\\<e^ . (1.31) 



Alternatively, in analogy to eq. ( |1.17| ), we may group the fault paths depending on the earliest faulty 
operation; 

L 

l^^faulty^^l^^^^)^^ (1.32) 

where 1 1 \d(JJr))\ | < 5. We conclude that 

6 <2Le . (1.33) 
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so that we obtain for local non-Markovian noise a similar result as for local Markovian noise — 
in both cases, the accuracy 1 — S of the noisy quantum circuit decreases at most linearly with the 
circuit size L. This illustrates that fully coherent noise, where the environment can store and process 
quantum information allowing different fault paths to interfere quantum mechanically, does not 
alter our view of quantum computation as a discrete model of computation similar to the model 
of modern digital computers (or, more abstractly, classical Turing machines). As we shall discuss 
in the following sections, methods of quantum fault tolerance can replace e in eq. ( |1.33| ) by an 



arbitrarily small effective noise strength showing that, just like Markovian noise, non-Markovian 
noise is not in principle an obstacle to large-scale quantum computation. 

1.2.3.1 Assessment and examples 



We note that eq. ( |1.33| ) was derived based on two essential assumptions: First, we assumed that 
the final noisy quantum state can be expanded perturbatively as a sum over fault paths. Secondly, 
we assumed that fault paths with many faults are suppressed in the sense of eq. ( |1.31| ). We can in 
fact define local non-Markovian noise in terms of these two assumptions, even if the noisy unitary 
evolution is not strictly local. We then generally say that 

Definition (Local non-Markovian noise). Noise is local and non-Markovian if the noisy evo- 
lution can be expanded as a sum over fault paths, where faults are described as (differences of) 
unitaries acting between the quantum computer and the environment and the norm of the sum of all 
the fault paths with faults in any r specific elementary operations is upper bounded by e^ for some 
constant noise strength s. 

This definition is very similar to our definition of local Markovian noise in the previous section: 
In both cases, the noisy evolution is expanded as a sum over fault paths and also noise is weak in the 
sense that, as the total number of faults in a fault path increases, the fault path norm is suppressed 
exponentially. In addition, in both cases noise can be correlated both in space and in time since 
we place no restrictions on the form of the fault operators that appear in each fault path, which are 
allowed to be arbitrarily and even adversarially correlated. 

The important distinction between the two cases is that while, for local Markovian noise, differ- 
ent fault paths do not interfere, for local non-Markovian noise, we make the worst-case assumption 
that all the fault paths do interfere coherently and the environment is not traced over until the end of 
the quantum computation. Therefore, while for local Markovian noise the strength e can be viewed 
as a probability (the probability for the occurrence of a single fault), for local non-Markovian noise 
the strength e corresponds in essence to a quantum amplitude (the amplitude for a term with a single 
fault in the final quantum state). If we were to use distinct symbols for the two cases, s and e' respec- 
tively for Markovian and non-Markovian noise then, since probabilities are squares of amplitudes, 
we expect e ^ {e') ; thus, requiring that the fault amplitude e' is small (say, less than 1.0 x 10~^) 
implies that the fault probability e is even smaller (in this case less than 1.0 x 10~^). 

The local non-Markovian noise model describes several noise processes for which the environ- 
ment interacts coherently with the quantum computer over long time scales; below, we discuss four 
examples of noise processes that give rise to local non-Markovian noise. 



Local Hamiltonian noise. The ideal unitary evolution in eq. ( |1.25| ) is generated by a time-dependent 
Hamiltonian 

Hqc = 2_\ -^QC ' ^^^^ that Uj = exp f —UoHq^ ) . (1.34) 

j 

The Hamiltonian that describes the noisy evolution of the quantum computer and the environment 
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then has the general form 

H = Hqc + He + Hqce , (1.35) 

where He generates the evolution of the environmental degrees of freedom, and Hqce describes 
the interaction of the quantum computer and the environment that introduces noise. If noise is local, 
Hqce has the same locality as Hqc, i.e., 

HQCE^Y^^hcE^ (1-36) 

3 

where the support of Hq^^ is contained in the union of the support of Hq^ and the environment. 
We can study the noisy evolution generated by H during a time interval of duration to pertur- 
batively; if we divide this interval into N micro-intervals each of duration Ato = t^/N, then 



N 



UQCEito, 0) = lim n U^cU^ceU^ , (1-37) 



n=l 



where Uqq, Uq^j^, and U^ denote the evolution during the n-th micro-interval according to Hqc, 
Hqce, and He respectively. After expanding^ 



U^CE^Ui^-'^^o^CE) ' (1-38) 

j 

and substituting in eq. ( |1.37| ), we can open the parentheses to obtain a perturbative fault-path expan- 
sion. The noisy implementation of the ideal unitary Uj then takes the form 

Uf'^ = {lQC + Fj)Uj; (1.39) 

here, Iqc denotes a sum of all the fault paths where, in every micro-interval, we insert either the 
identity or a micro-fault fj> acting on an operation with label j'^ j where 

/,v = -iAtoH^QCE ■ (1-40) 

It follows that Fj includes all remaining fault paths where a micro-fault fj acting on the operation 
with label j is inserted in at least one micro-interval. 

We can express Fj as a sum of terms labeled by the micro-interval where the earliest micro-fault 
fj is applied on the operation with label j; if we denote by llQCEi^tQ, Atg) the entire evolution 
generated by Hqce between the i-th and j-th micro-intervals with i < j, then 

^^^i^^T. f ^Qc^(^^^'^^s^')n ( ^gc/f'^ n (^+/.o^i 1 1 u^j • (1-41) 



Now, since the operator norm is unitarily invariant, each term in the sum in eq. ( |1.41[ ) has norm 

WfjWoo = Ato\\H^QCE\\oo, ^nd so \\Fj\\oo < to\\H^QCE\\oo' 

In fact, we can perform a similar perturbative expansion to analyze faults in any specific subset 
of the L elementary operations in the quantum circuit. It follows that local Hamiltonian noise 
satisfies eq. ( |1.31| ) with 



e = to -maxlliJ^^^lloo • (1.42) 

As expected, the noise strength e depends on the strength of the interaction term Hq^^ between 
the quantum computer and the environment and also the time during which is interaction is acting. 
On the contrary, we observe that e does not depend on the strength of the term He which describes 

^ Since we take the limit A/"— ^oo, we have Ato— ^0 and we may keep only the linear term in the expansion. 
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the internal evolution of the environment. Moreover, we note that while Hq^^ is assumed to act 
locally on the quantum computer in the sense of eq. ( |1.36| ), the derivation of eq. ( |1.42| ) did not rely 
on making any assumptions about He which is completely arbitrary. 

Long-range static noise. In certain systems, noise can arise due to static — i.e., time-independent — 
interactions among pairs of qubits of the quantum computer, where these interactions do not depend 
on the ideal circuit that is being implemented. Such non-local noise can be modeled by the Hamil- 



tonian in eq. (1.35) with 



where i?(j,it) is supported on qubits j, k and the environment and we sum all unordered pairs (j, k) 
of qubits. 

We can perform a similar perturbative expansion as in the case of local Hamiltonian noise, 
except that now the two qubits in the support of any micro-fault may not be directly interacting via a 
unitary gate. Despite this difference which necessitates a more complicated combinatorial analysis 
(see the references), it can be shown that long-range static noise satisfies eq. \\3\) with 



5= /cto max^||i7(_^-^;,)||oo , (1.44) 

where c = 2e provided e^ < e and it is understood that, if Hqce is time dependent, the maximum 
is also taken over all times. 



Gaussian noise. In a variety of physical setting where the qubits of the quantum computer are 
coupled to a large number of environmental degrees of freedom, the environment can be well 
approximated as a collection of uncoupled harmonic oscillators obeying Gaussian statistics; the 
Hamiltonian of the environment is 

He ^^oJkdlak , (1.45) 

k 

where ak are bosonic annihilation operators satisfying [a/e,a[/] = Skk'- In this spin-boson model 
of the noise, the interaction between the quantum computer and the environment is described by a 
coupling of each qubit to a linear combination of oscillator amplitudes 

Hqce = ^ crm{x) ^m{x, t) (1.46) 

with 

c^m{x, t) = e^^^^0^(x, t)e-^^^^ = Y. [dkA^. t)ake-''^^ + gl^{x, t)4e^^"^) , (1.47) 

k 

where x labels a qubit's position and a^(x) with m G {x, y, z} are the three Pauli operators on the 
qubit with label x. 

The statistics of the environment amplitudes (l)m{x^ t) are Gaussian in the sense that the n-point 
correlation functions with respect to the environment state |0o)^ vanish for n odd, while for n even 
they obey Wick's theorem: 

{(t^mi{xi,ti) • • • (l)mr,{Xn,tn)) = ^ A(n, Z2) ' ' ' A(z^_l, Z^) , (1.48) 

(^1,^2), ••• ,(^n-l,^n) 

where A(p, q) = {(j)mp{xp^ tp)(j)rn^{xq^ tq)), and we sum all ways of dividing the label 1 to n into 
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n/2 unordered pairs. By performing a perturbative analysis similar to the case of long-range static 
noise (see the references), it can be shown that Gaussian non-Markovian noise satisfies eq. \\3\) 
with 



e= cm^^ f f J2 l^(l'2)|, (1.49) 

where the first integral denotes an integration over the qubits in the support of the unitary Uj and the 
time interval during which this gate is implemented, and the second integral denotes an integration 
over all the qubits of the quantum computer and the total duration of the quantum computation. 

Local leakage noise. The qubits of the quantum computer are in practice always realized as two- 
dimensional subspaces inside a multi-dimensional system; the Hilbert space Hqc of the quantum 
computer then has a natural extension to 

n'^c = 'HQC(B'H^c, (1-50) 

where the leakage space ^^qc includes all states outside the two-dimensional qubit subspaces — in 
most settings, 1-Lq(j is a tensor product over leakage spaces corresponding to each qubit. 

Now, the Hamiltonian that describes the noisy evolution of the quantum computer and the en- 



vironment has the same general form in eq.(1.35). 



H = Hqc + Hqc + He + Hqce , (1.51) 

where Hq^j generates the evolution in the leakage space, and Hq^^ describes the interaction be- 
tween the extended space of the quantum computer (the qubits and their leakage spaces) and the 
environment. If noise is local, Hq^^ has the same locality as Hqc, i.e., 

Hqce = /Jh^ce 5 (1-52) 

3 

where the support of ^Hqq^ is contained in the union of the support of Hq^, the leakage space, 
and the environment. By repeating the same analysis as for local Hamiltonian noise, we find that 
local leakage noise satisfies eq. \\3\) with 

5 = to-max|pi/§^(^^||oo. (1.53) 

3 



1.3 Encoded quantum computation 

When we desire to implement long computations — i.e., when the size L of the quantum circuit is 



large — an accuracy that decreases linearly with L as in eqs. \\.\6) and ( |1.33| ) is not satisfactory. Of 
course, we could achieve an accuracy independent of L if s were a decreasing function of L thus 
making Le a constant, but this is certainly not a physically reasonable assumption — we cannot hope 
that the hardware of the quantum computer will get less and less noisy the longer we keep quantum 
computing! 

In order to obtain the results of a quantum computation with constant accuracy, some method 
for detecting and correcting the errors that are introduced by the noisy hardware is necessary; we 
say that such a method of computation is fault tolerant. The basic idea of fault-tolerant computation, 
whether classical or quantum, is the use of redundancy. Every hardware operation in the circuit to 
be implemented is replaced by a collection of several hardware operations which are designed to be 
more robust to local noise than a single hardware operation alone. 
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Fig. 1.3. On the left, a NOT gate is applied on the bit value carried in a single wire — e.g., if the NOT operation 
is implemented by a CMOS inverter, the input wire controls the gate voltage and the output is taken as the 
drain voltage. The input and output wires on the left are replaced on the right by cables each comprising n 
wires, and also the NOT gate on the left is replaced on the right by an organ comprising n parallel NOT gates 
(here, n = 3). 

A formal method for introducing redundancy is via the use of error-correcting codes. For clas- 
sical computation, the simplest example of a redundant encoding of information is based on the 
repetition code: "To protect information traveling from gate to gate, we replace each wire of the 
noiseless circuit by a cable of n wires (where n is chosen appropriately); each wire within the cable 
is supposed to carry the same bit of information, and we hope that a majority will carry this bit even 
if some of the the wires fail.'Q To protect information during the execution of each gate, we also 
replace each gate in the noiseless circuit by an organ comprising several gates. The organ operates 
on the information carried by the wires inside the cables in the same way that the initial unencoded 
gate operated in the information carried by single wires — e.g., a NOT gate must be replaced by n 
NOT gates acting in parallel on every wire in a cable as in fig. |1.3[ and similarly a AND gate must be 
replaced by n AND gates acting in parallel on corresponding pairs of wires in two cables. Organs 
also include a procedure for detecting and correcting faults in the noisy hardware; as discussed by 
von Neumann (see the references), this may be implemented for each wire by copying the value of 
every bit to a larger number of k bits, randomly permuting all the resulting n-k bits, computing the 
majority function in parallel n times on disjoint sets of k of these permuted bits, and having the n 
outputs form the output wire. 

Similar redundant encodings are also possible for quantum information where, for historical 
reasons, cables are now called blocks and organs are called gadgets^ Each qubit in the noiseless 
quantum circuit is replaced by a block of n encoded qubits; the joint state of the encoded qubits is 
supposed to carry the same quantum information as the state of the initial qubit, and we hope that 
this quantum information can be recovered even if faults occur on some of the encoded qubits. In 
addition, each elementary operation in the noiseless quantum circuit is replaced by a gadget that 
comprises several elementary operations acting on the encoded qubits in a block or across multiple 
blocks. A gadget is designed to operate on the quantum information carried by the encoded qubits 
in the same way that the initial unencoded operation acted on the quantum information carried 
by single qubits, and it also includes a procedure for detecting and correcting faults in the noisy 
hardware. The encoding of quantum circuits is shown schematically in fig. |1.4| 

But what do we mean when we say that the joint state of the encoded qubits carries the same 
quantum information as the state of the initial unencoded qubit? And how do we hope to recover the 
quantum information that is carried in a block if local noise acts on the encoded qubits? To answer 
the first question, consider a pure state \iP)biR supported on Hbi ® T-Lr where Hbi is the Hilbert 
space of a qubit Bi and T-Lr is the Hilbert space of a reference system R. The encoded version of 
|^)Bii? is a pure state \iI^)br supported on %^ Hr, where %^ is the code space of a block B of 
n encoded qubits Si, ... , Bn. We say that the states \iP)biR and \iP)br carry the same quantum 
information because \iP)br can be obtained from \iP)biR by applying an isometry that maps Hbi 
to T-L^; or in other words, there exists a unitary decoding unitary operator t/dec acting on the block 

^ Quote from Gacs's Reliable computation; see the references. 

^ People also use the term rectangle instead of gadget when they think of the latter's pictorial representation as in Fig. 1.4 
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time 



Fig. 1 .4. Two elementary operations (single-qubit preparations, unitary gates, or measurements) in the noise- 
less quantum circuit on the left are replaced in the actual noisy quantum circuit on the right by two gadgets. 
For illustration, the elementary operations comprising each gadget are enclosed by dashed-line rectangles and 
only two elementary operations per gadget are shown. 



such that 



Tr 



{RBi)- 



-{Uaec^InW 



BR 



/BiR 



(1.54) 



where Ir is the identity operator on the reference system and Tr/^^^)± denotes a trace over every- 
thing else except for the reference system and the qubit Bi. 

To answer the second question, we recall from Section [TT] that the basic idea of quantum error 
correction is error discretization. To monitor the effects of noise, we partition the entire Hilbert 
space of the encoded qubits into mutually orthogonal subspaces and, if there is no noise, we de- 
mand that the support of \tp)BR coincides at all times with the code space T-L^ (and the reference 
system). In the presence of noise, our strategy is to detect periodically whether \^P)br develops a 
non-zero overlap with any other subspace H^ labeled by a non-trivial syndrome s (cf., eq. ( |1.2[ )), 
in which case we apply a recovery operation that returns the support of |^)^i^ to T-L^. Physically, 
distinguishing on which subspace the state \iIj)br is supported can be implemented by performing 
a generalized measurement A^{^} jointly on the encoded qubits, processing the measurement out- 
come 11 to determine the syndrome s and hence the subspace H^ on which \iP)br has been projected 
by the measurement^ and applying an operation TZg on the encoded qubits conditioned on the value 
of s that maps H^ to the code space H^. 



We may revise eq. ( |1.54| ) to include cases when the quantum information has been afflicted by 
noise: We now say that \iP)biR and the noisy \iP)br carry the same quantum information if 
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^Ir , 



(1.55) 



(1.56) 



where Udec^^R are the physical operations corresponding to applying the unitaries Udeci Ir^ respec- 
tively, and we have suppressed the classical on-the-side computation which determines s given /x. 
It is worth of notice that the classical bits carrying the outcome ii of the syndrome measurement. 



which are traced over in eq. ( |1.55| ), carry information about the subspace on which the quantum in- 
formation was encoded prior to the decoding; for this reason, they are often referred to as syndrome 
bits. 

The combined operation V is called a decoder. Because decoders output qubits which are un- 
encoded and therefore unprotected from noise, we will never use a decoder in our actual noisy 
quantum circuits — encoded quantum information will never be decoded. A noiseless ideal decoder 
is however very useful as a tool for formalizing the requirement that a gadget operates on the en- 



^ This processing may be performed in a classical on-the-side computer. 



18 



rR 



M^)Blfi 



IBR 



B2 


^XK. 




^MV. 




wv 


B3 


•^IjuI 




v.. 




Z^dec 


Bn' 

V 


XH'f ' -J 

Km K^J 




^XK. 



V 



\ 



^i 



m 



Fig. 1.5. A decoder V comprises a syndrome measurement M{^y, followed by a classical computation 
that outputs the syndrome s with input the measurement outcome /i, followed by a recovery operation 1Zs 
conditioned on the syndrome, followed by a decoding operation U^ec mapping the code space of the block of 
encoded qubits 5i , ^2 , • • • , ^n to the Hilbert space of the qubit Bi . Our notation for V shows that the input 
consists of n encoded qubits and the output of one decoded qubit, qubit Bi . As denoted by the bold gray 
line, the decoder output also includes the qubits ^2 , • • • , ^n and the syndrome bits that carry the outcome of 
the syndrome measurement; the state of these systems ends up being always a tensor product with the joint 
state of qubit Bi and the reference R. When we speak of an ideal decoder, which is an imaginary noiseless 
operation, we draw V and the operations comprising it inside wavy boxes. 



coded qubits in the same way that the operation that was replaced by the gadget acted on the initial 
unencoded qubits. 

It is convenient to denote noiseless ideal operations using wavy boxes and noisy operations 



using square boxes; e.g., an ideal decoder is shown in fig. |1.5[ With this notation, (a) an operation 
V\^) that prepares the single-qubit pure state |^) in the noiseless quantum circuit is replaced in the 
noisy quantum circuit by a gadget implementing the operation VP,s^ on the encoded qubits where 
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(1.57) 



(b) an operation U that applies the single-qubit unitary U is replaced by a gadget implementing the 
operation Ul on the encoded qubits where 



Ul 



(\y\Ji 



V 



r\>^ 



\ 



V 



U 



(1.58) 




(and similarly for multi-qubit unitary operations), and (c) an operation M.\a] that measures a single- 
qubit observable a with eigenvalues {a} is replaced by a gadget implementing the operation A^ j^j 
on the encoded qubits where 
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(1.59) 



Eqs. ( |1.57| ), ( |1.58| ), and ( |1.59| ) capture what we mean when we say that a gadget implementing the 
operation Ol simulates an unencoded operation O. 

But apart from simulating the desired operation on the encoded qubits, gadgets need also con- 
tain a mechanism for detecting and correcting faults that may occur in the elementary operations 
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comprising them. We can introduce such a mechanism by inserting an error recovery operation 
inside every gadget. Specifically, the gadget that implements Ul now becomes 



Ul 
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Ul 
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STl 



(1.60) 



where ETZ comprises a syndrome measurement followed by a recovery operation, 
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(1.61) 



w ' 



For the gadget that implements Vhy we insert similarly an error recovery operation following Vhy 
while the gadget that implements M^^^ is not modified since its output is a classical number, the 
measured eigenvalue a, which is protected from noise provided the classical computer that processes 
the measurement outcome is robust. 

We do not have enough time and space to discuss explicit gadget constructions here. What 
is important at this abstract level is that each gadget simulates an operation in the noiseless quan- 
tum circuit, and that we hope that gadgets can be more robust to noise than unencoded operations 
because faults inside each gadget can be detected and corrected by the error recovery operation. 



13,1 Properties of noisy gadgets 

To make progress, we need a method for formalizing the degree to which gadgets are protected from 
noise. We recall from Section [T.2.1| that if the noisy quantum circuit is afflicted by local noise, we 
can expand the noisy evolution perturbatively as a sum of fault paths; each fault path identifies a 
specific subset of all the elementary operations where faults have been inserted, while there are no 
fault insertions in all remaining elementary operations. The idea of fault-tolerant constructions is to 
ensure that gadgets operate reliably for all the fault paths with no more than a certain number t>l 
of insertions of faults inside them; the intuition is that each gadget will then be more robust to noise 
than any single elementary operation alone because the first contribution to a gadget's failure comes 
at order t+1 of our perturbative fault-path expansion. 

1.3.1.1 Good gadgets 

But what do we mean when we say that gadgets operate reliably for fault paths with at most t faults 
inside them? To formalize this requirement we need to consider each gadget together with the error 
recovery operations of the immediately preceding gadgets; we will refer to a gadget together with 
its preceding error recovery operations as an extended gadget. The idea is to construct gadgets such 
that (a) for each fault path with at most t insertions of faults inside a measurement extended gadget, 
the noisy gadget is equivalent to an ideal measurement of the gadget's ideally decoded input: 
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(1.62) 



(b) for each fault path with at most t insertions of faults inside an extended gadget simulating a 
single-qubit unitary gate, the gadget is equivalent to applying the ideal unitary gate to the gadget's 
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(1.63) 



<t insertions 

and (c) for each fault path with at most i insertions of faults inside a preparation gadgej^ the gadget 
is equivalent to the ideal preparation: 
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(1.64) 



Here, we have illustrated insertions of faults as couplings between the noisy elementary operations 
inside the gadgets and the environment denoted by a think red line. 



Figuratively, eq. (1.62) allows us to create ideal decoders out of measurement extended gadgets 



which contain at most t faults, eq. (1.63) allows us to propagate ideal decoders to the left through 



unitary-gate extended gadgets which contain at most t faults, and eq. ( |1.64| ) allows us to annihilate 
ideal decoders in preparation gadgets which contain at most t faults. As ideal decoders are created 
in measurement gadgets, propagated through unitary-gate gadgets, and annihilated in preparation 
gadgets, the noisy encoded quantum circuit is transformed to a noiseless unencoded quantum circuit. 
We conclude that for all the fault paths with at most t faults in each and every gadget, a noisy 
quantum computer is equivalent to a noiseless ideal quantum computer, in the sense that they both 
produce the same probability distribution for the final measurements which determine the computa- 
tion result. However, for fault paths for which there are more than t faults inside one of the extended 
gadgets, the output probability distributions from the noisy and the noiseless quantum computer are 
not guaranteed to agree; thus it appears that the accuracy 1 — 5 of a noisy encoded quantum circuit 
composed of fault-tolerant gadgets scales with 



5 ^ 5^+1 , 



(1.65) 



which should be compared with the scaling 5^5 for a noisy unencoded quantum circuit — cf., 
eq. \\.\6\ for local Markovian noise and eq. ( |1.33| ) for local non-Markovian noise. To prove that 



this scaling indeed holds and to determine the proportionality constant in eq. ( |1.65| ) requires that we 
analyze what happens when a gadget is afflicted by more than t faults. 

1.3.1.2 Bad gadgets 



Eqs. ( |1.62| ), ( |1.63| ), and ( |1.64[ ) show that extended gadgets that contain at most t faults can be viewed 
as implementing ideal operations acting on unencoded qubits. But if more than t fault occur inside 
a gadget, there are no guarantees what may happen. Although we cannot in general say much about 
what actually happens unless we know details about the noise and how gadgets are constructed, 
we will be satisfied if we can show that extended gadgets which contain more than t faults can be 
viewed as implementing some noisy operations acting on unencoded qubits. 

The basic tool we will need is a decomposition of the identity, i.e. do nothing, operation in 
terms of an ideal decoder-encoder pair : 

^ And similarly for multi-qubit unitary gates where we include all preceding error recovery operations. 

^ The extended preparation gadget coincides with the preparation gadget since preparation gadgets have no input and, therefore, they 
are not preceded by any error recovery operations. 
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Fig. 1.6. A reversible ideal decoder, to be contrasted with the irreversible ideal decoder in fig. |1.5| First, 
the syndrome measurement A^{/i} is replaced by a unitary operation Udet that acts on a larger Hilbert space 
which includes ancillary syndrome qubits; this unitary maps coherently the state of the syndrome qubits to a 
state corresponding to the subspace s on which the input block is supported. Secondly, the recovery operation 
7ls is replaced by a controlled unitary operation Urec that applies the appropriate recovery unitary on the block 
depending on the syndrome s carried by the syndrome qubits. In the figure, we assume there are m syndrome 
qubits and that each one is prepared in the state |0). 




(1.66) 



where the operation V ^ is an ideal encoder. Of course, for ideal decoders to have an inverse, they 
need to be implemented by a reversible circuit. In fact, we can easily modify our definition of ideal 



decoders in Section [L3] to achieve reversibility. Fig. |1.6| shows such a reversible ideal decoder; an 
ideal encoder can then be implemented by simply executing the circuit in this figure backward in 
time. 

Decoder-encoder pairs are trivial (they multiply to the identity operation) but they can be useful 
for understanding the properties of our noisy circuits if we insert them strategically at appropriate 
places. As a start, consider how to modify eq. ( |1.63| ) in the case there are more than t faults. The 
following property now holds: For each fault path with more than t insertions of faults inside an 
extended gadget simulating a single-qubit unitary gate, the gadget is equivalent to applying some 
noisy gate to the extended gadget's ideally decoded input: 
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where the noisy gate is 
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There are two points about eq. \\.61) deserving emphasis. First, this property replaces the 



entire extended gadget by a (noisy) gate rather than replacing the gadget alone; cf. eq. (1.63) and 



22 



see Section [T. 3. 1.3 1 below. Secondly, the noisy gate that replaces the extended gadget depends on 
both the faults inside the extended gadget and also on the syndrome bits (found inside the bold gray 



lines) which are input to the ideal encoder. As our notation in eq. ( |1.68| ) is intended to illustrate, 
the syndrome bits can be viewed as a fictitious environment which operates together with the actual 
environment associated with the noise. 

Similar properties as eq. ( |1.67| ) can be derived for measurement extended gadgets and prepara- 
tion gadgets; schematically, 
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where the noisy measurement on the right-hand side can be obtained by inserting a decoder-encoder 
pair on the left-hand side, and 
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1.3.1.3 Truncation 



Figuratively, eqs. ( |1.67| ), ( |1.69| ), and ( |1.70| ) allow us to create, to propagate to earlier times, and to 
annihilate ideal decoders in the case when extended gadgets are bad containing more than t faults. 



These properties are therefore complementary to eqs. ( |1.62| ), ( |1.63| ), and ( |1.64| ) that apply to good 
extended gadgets containing at most t faults. 

Of course, one difference between the two cases is that in one case noisy unencoded operations 
appear on the right-hand side, while in the other case the unencoded operations are noiseless and 
ideal. But here we would like to discuss a second difference that was mentioned already in the 
previous section; namely, while eqs. ( |1.67| ), ( |1.69| ), and ( |1.70| ) replace the entire bad extended gad- 
gets by some (noisy) unencoded operations, eqs. ( |1.62| ), ( |1.63| ), and ( |1.64[ ) only replace the gadgets 
contained inside the good extended gadgets by the ideal unencoded operations. This modification 
is important in order to prevent overcounting faults in successive, and therefore overlapping, bad 
extended gadgets. 

To be concrete, imagine that we are to encode a quantum circuit comprising just a single-qubit 
preparation V\^) followed by a single-qubit measurement M^a}- The encoded quantum circuit 
comprises two overlapping extended gadgets: 



preparation 

extended 

gadget 






sn 



^w 



measurement 

extended 

gadget 



(1.71) 



where the error recovery step is contained in both extended gadgets enclosed in the square gray 
boxes. Now, consider a fault path with more than t insertions of faults in each extended gadget. 
Considering the two extended gadgets in isolation, we may think that the noisy encoded circuit 
equals some noisy preparation followed by some noisy measurement: 
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However, we soon realize that this may not always be a satisfying answer. Because the two extended 
gadgets overlap, it is possible that the total number Nf of faults in the given fault path is less than 
2(t+l); thus the two noisy unencoded operations on the right-hand side appear at order e^f in 
our perturbative fault path expansion, which is less than the order (5^+^) we would expect based 
on the fact that each extended gadget by itself fails at order e^^^. Clearly, the problem is that our 
naive estimate (5^+^) double counts each fault inside the error recovery step shared by the two 
overlapping extended gadgets. 

Nevertheless, this complication is a red herring, and our naive estimate can actually be jus- 
tified: To formally obtain eq. ( |1.72[ ), we first need to use eq. ( |1.69| ) thereby replacing the entire 
measurement extended rectangle by a noisy unencoded measurement: 
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>t insertions 

where we observe that the error recovery step has been removed from the preparation extended 
gadget — we say that the gadget has been truncated. To annihilate the ideal decoder at the next step, 
we need to consider how many faults are contained inside the truncated preparation gadget; if there 



are at most t faults then we can use eq. (1.64) to obtain 
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whereas, if there are more than t faults, we can use eq. ( |1.70| ) to obtain eq. ( |1.72| ); see, however, 
the next paragraph which explains why eq. ( |1.70| ) applies to truncated extended gadgets. Since in 
this second step what matters is the number of faults in the truncated preparation gadget, there is 
no double counting of faults; thus the two noisy unencoded operations on the right-hand side of 
eq. (1.72) do appear at order (5^+^) in our perturbative fault path expansion, as desired. 

We need to note that eqs. ( |1.62| ), ( |1.63| ), and ( |1.64| ) for good extended gadgets and eqs. ( |1.67| ), 
( |1.69| ), and ( |1.70| ) for bad extended gadgets, which have been formulated for the full extended 
gadgets, apply in the same way to truncated extended gadgets. Indeed, the ideal decoders contain a 



noiseless ideal error recovery step — cf . fig. ( |1.5| ) — which can be used to replace the truncated noisy 
error recovery steps, thereby reassembling the full gadgets for which the properties apply. 

Although we have illustrated the concept of truncation with a simple example, a similar trunca- 
tion procedure can be used for any fault path which leads to more than two successive bad extended 
gadgets, including gadgets that simulate unencoded operations on more than one qubit: Starting 
from the latest bad extended gadgets which are not succeeded by any other bad extended gadgets, 
we progressively move to all earlier bad extended gadgets, one gadget at a time. At each step, we 
truncate the bad extended gadget under consideration from the error recovery steps it shares with all 
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its succeeding bad extended gadgets (which may themselves be truncated or not), and we label the 
truncated gadget as good or bad depending on the number of faults it contains after the truncation 
(if it contains at most t faults, it is declared good, otherwise it is declared bad). Eventually we 
reach the earliest bad extended gadgets, and we truncate the good extended gadgets that preceded 
thenq^ In the end, the successive bad extended gadgets are divided into non-overlapping truncated 
extended gadgets which have been declared good or bad depending on the number of faults they 
contain and which can be replaced by either noisy or ideal unencoded operations respectively by 
using the corresponding gadget properties. 

1.4 Coarse-grained noise and level reduction 

We have now assembled all the properties we need to characterize noisy gadgets. If we combine 
all the pieces together, we can arrive at a helpful description of the noise acting on the encoded 
quantum computer. 



We recall from Section 1.2.1 [that, without encoding, the only fault path leading to the noiseless 



ideal evolution is the trivial fault path which contains absolutely no faults; we are then forced to 
decompose the noisy evolution into an ideal and a faulty part as in eq. ( |1.11| ) for local Markovian 



noise and eq. ( |1.28| ) for local non-Markovian noise. In contrast, in an encoded quantum computa- 
tion, many more fault paths lead to the noiseless ideal evolution; now, for local Markovian noise, 
we can decompose the noisy evolution as 

^noisy _ ^good ^ ^bad ^ (1.75) 

and similarly for local non-Markovian noise, 

where the unnormalized density matrix pS^^^ and the unnormalized pure state |^)s°°*^ are sums of 
all the fault paths with at most t faults in each and every extended gadget, while (^^^ and \^)^^^ 
sum of all the remaining fault paths. 

In what sense are the fault paths included in (^^^^^ and |i?)g^^^ good? For each (good) fault 
path with at most t faults in each and every extended gadget, eqs. ( |1.62| ), ( |1.63| ), and ( |1.64| ) apply. 



We now consider the entire encoded quantum circuit and, by using these properties, we first create 
ideal decoders in the measurement gadgets. Then, we propagate the ideal decoders to earlier times 
through unitary-gate gadgets. Finally, we annihilate the ideal decoders in preparation gadgets. As 
the ideal decoders appear, move to earlier times, and finally disappear, the entire encoded quantum 
circuit afflicted by the given good fault path is shown to be formally equal to the ideal quantum 
circuit that the gadgets simulate. Since this is true for every good fault path separately, by linearity 
it is also true for the sum of all of them; thus (^g^^^ and \i})^^^^ lead (after normalization) to the 
same probability distribution for the final computation result as a noiseless ideal quantum computer 
would. 

We can now estimate the accuracy 1 — 6 of the encoded quantum computation. For local Marko- 
vian noise, 6 can be bounded as in eq. ( 1.13| ) by the norm of the difference of the final noisy super- 



operator minus its (normalized) good part: 

^good / 1 \ 

' S IIP"'- - ^f^ll. i (l + T^W^J IIC"II. . 0.77) 

and we have used the triangle inequality multiple timeq^ Similarly, for local non-Markovian noise, 

^ These extended gadgets remain good after the truncation as they were aheady good prior to it. 

2 We have (5 < ||apSood ^ ^bad||^ < |a| + ||C^^d||i, where |a| =Tr-i(pg°°d) -1 < ||C^^^||i/(l - IK^^'^Hi). 
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5 can be bounded as in eq. (1.29) by the norm of the difference of the final noisy pure quantum state 



minus its (normalized) good part: 

* s mtr-' - ^p^ii < 2 (i + T^^) iii*>'-^ii . (1-78) 

It remains to obtain upper bounds on (^^^"^ and 1 19) ^^*^ which are sums of all {bad) fault paths with 
more than t faults in at least one gadget. For each bad fault path, the noisy encoded quantum circuit 
can be analyzed by using the gadget properties: We first consider whether each extended gadget 
contains at most t faults or more than t faults, declaring the former gadgets good and the latter bad. 
If multiple successive extended gadgets are declared bad, we use the truncation procedure described 
in Section ( |1.3.1.3| ) to divide them into non-overlapping truncated extended gadgets which are good 



or bad depending on the number of faults they contain. Eventually, we use eqs. (1.67), (1.69), and 



( |1.70| ) for the good extended gadgets (truncated or not) and eqs. ( |1.62| ), ( |1.63| ), and ( |1.64[ ) for the bad 



extended gadgets (also, truncated or not). Every good extended gadget is thereby replaced by the 
noiseless ideal operation the gadget simulates, every bad extended gadget is replaced by some noisy 
operation, and thus the noisy encoded quantum circuit as a whole is replaced by a noisy unencoded 
quantum circuit. 

If we now let C^^^ denote the set of all L extended gadgets in the encoded quantum circuit, then 
by analogy to eq. ( |1.14| ) we may write 

L 

c'^' = E(-ir"' E c(c«), (1.79) 

where the second sum is over all subsets Ci- of C^^^ of cardinality r, and C,{Cr ) denotes a sum 
of all the fault paths for which all the extended gadgets in Cr are declared baq^ 

To gain intuition about how to proceed, consider the simplest case r = 1 when Q{C\ ^) is a sum 
of all the fault paths for which the single extended gadget in a specific set C\ ^ is bad. Since for an 
extended gadget (truncated or not) to be bad it needs to contain more than t faults, we can generalize 



eq. ( |1.14| ) to obtain 



where Lq is the number of elementary operations in the extended gadget in C} \ the second sum 
is over all subsets Cg of s of these operations, and C(C's) is a sum of all the fault paths with faults 
applied on all operations in Cg [5 

For the general case of r bad extended gadgets, it suffices to perform a similar inclusion- 
exclusion analysis independently in each gadget: 



ac^'')-\{\j:i-^y^-'-'i^^^^^ E c(c.,) 



(1.81) 



where Lq now denotes the number of elementary operations in the largest extended gadgejj Cr ( j) 

^ Whether each extended gadget in Cr is truncated depends on the fault path; however, for each specific fault path, we can first 
use the truncation procedure to decide which extended gadgets need to be truncated, and then we can unambiguously declare every 
extended gadget (truncated or not) as being either good or bad. 

^ We first sum all C,{Ct+i) accounting correctly for all the fault paths with exactly t-\-l faults; however, all the fault paths with exactly 
t-\-2 faults are overcounted (^i-,^) — 1 = {^\ ) times. So next, we subtract the sum of all (^(Ct+2) multiplied by (*"^ ) , but in doing 
so we undercount all the fault paths with exactly t+3 faults (*+J) - (*+-^) (*t2) -1 = - {^^\~^) times. And so on. 

'^ The number of elementary operations may vary among gadgets. In addition, some gadgets may be truncated depending on the fault 
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Fig. 1.7. The encoded operations in an encoded quantum circuit are executed by using gadgets (here, one 
gadget is shown, along with two of the elementary operations it contains); the elementary operations inside 
each gadget are afflicted by local noise (either Markovian or non-Markovian) of strength e. The physical 
noise can be coarse-grained into a local effective noise with renormalized strength e*^^^ of order 6^+^ acting 
on the gadgets themselves. The noise coarse-graining allows us to concentrate on the gadgets alone and forget 
about the elementary operations inside them — we say that the encoded quantum circuit is level reduced to an 
equivalent unencoded quantum circuit, where the effect of coding is to map the physical noise strength e to 
the effective noise strength e^^\ 



denotes the set of elementary operations in the j-th bad extended gadget in Cr , and Cs^ denotes a 
subset of Sj of the elementary operations in Cr (j). 

But I \C{Csj ) 1 1 1 < ^^^ by the definition of local Markovian noise, and we find 



Lq 



iic(c(^^)iii<n E (V ^^' 
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(1.82) 
(1.83) 



for some constant ^ > e(^o~^~^)^ (typically we are interested in small values e < l/(Lo— t— 1); 
then we may take ^ to be e). By replacing ( with |i9), we may repeat a similar calculation for local 
non-Markovian noise to obtain 



imc^'^m < [e^y 



(1.84) 



with £:(^) again as in eq. ( 1.83 ). 

Eqs. ( |1.82| ) and ( |1.84| ) tell us that if we choose any r extended gadgets, then the sum of all the 
fault paths for which all of them are bad has norm which is exponentially suppressed with r. This 
is exactly the condition we imposed for noise to be local, except this in this case we think of noise 
as afflicting the gadgets themselves instead of the elementary operations. While the strength of the 
noise acting on the elementary operations is e, the strength of the coarse-grained noise acting on the 
gadgets is £^^\ which scales as s^^^ because each gadget can tolerate up to t faults. We often refer 
to the encoded computation executed by the gadgets as a level-1 simulation of an unencoded level-0 
quantum circuit; in this language, what we have shown is that a noisy level-1 simulation afflicted 
by local noise with strength e can be viewed as a level reduced noisy level-0 simulation afflicted 
by a coarse-grained local noise with renormalized strength e^^^ (and the level reduction works the 



same for both Markovian and non-Markovian noise). Fig. |1.7| illustrates this coarse-graining level 
reduction procedure. 

path. By taking Lq to correspond to the largest extended gadget, we thus unavoidably include in the sum some extra fault paths 
which should not be counted. However, since eventually we will take the norm of both sides and use the triangle inequality, including 
these additional fault paths merely weakens our bounds. 
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The accuracy of an encoded computation afflicted by local Markovian noise can now be deter- 
mined from eq. \\.ll) by combining eqs. ( |1.79| ) and ( |1.82[ ); we find 



5 ^ llC'^^dlK < Y. (^) {^^"^y < {e-l)Le^^ , (1.85) 

r=l ^ ^ 

where, in the first step, we have kept only the leading order contribution. By comparing with 



eqs. ( 1.16 ) which corresponds to the case no encoding is used, we conclude that if e^^"^ < e then the 
encoding is in fact a good idea since it improves the accuracy 1 — 5 of the final computation result. 
The same conclusion also holds for local non-Markovian noise. 



1.5 The quantum accuracy threshold 

We are one breath away from the central result in the theory of quantum fault tolerance. You must 
have guessed the next step... If an encoded quantum circuit is more accurate than a quantum circuit 
which is not encoded, then why not apply the encoding to the encoded circuit itself, taking every 
elementary operation inside it and replacing it by a gadget; this doubly encoded quantum circuit 
should be even more accurate. In fact, why stop here? If we continue recursively re-encoding our 
encoded circuits, we expect their accuracy to steadily increase reaching any limit we please. 

To formalize this idea, let us consider the recursive construction of these multiply encoded 
quantum circuits we imagined above. At the base of our construction is the unencoded quantum 
circuit corresponding to our quantum algorithm; we say that this is our level-0 circuit in the sense 
that it does not use any coding. The next step is to replace every elementary operation in the level- 
circuit by the corresponding gadget; we say that this is our level- 1 circuit, performing a level- 1 
simulation of the level-0 circuit. Instead of physically implementing the level- 1 circuit as is, we may 
next replace every elementary operation in the level- 1 circuit by the corresponding gadget to obtain 



our level-2 circuit, and so on. Fig. [L8] illustrates this replacement procedure repeated k times; the 
final encoded quantum circuit, which is the one we do physically implement, performs a level-A: 
simulation of the level-0 circuit. 

The question is what is the accuracy of the level-A: circuit as a function of k. Estimating this 



accuracy is actually especially easy if we use the noise coarse-graining concept from Section [L4 
The noise afflicting the elementary operations in the level-A: circuit can be coarse-grained to give 
an effective noise that acts on the gadgets; if the physical noise is local and h as str ength e, the 
coarse-grained noise is also local and has renormalized strength e^^^ as in eq. (1.83). The noise 



coarse-graining level reduces the level-A: circuit to an equivalent level-(A:— 1) circuit; this level- 
(fc— 1) circuit produces the same probability distribution for the computation outcome as the initial 
level-A: circuit, but it is afflicted by an effective local noise of strength e^^\ Thus we have reduced 
the problem of estimating the accuracy of the level-A: circuit (the circuit we actually physically 
implement) to estimating the accuracy of the level-reduced level- ( A:— 1) circuit (which is, of course, 
imaginary as it represents the result of our noise coarse-graining procedure). 

We can next coarse grain the noise in the level- (A:— 1) circuit, thereby level reducing the initial 
level-A: circuit to a level-(A:— 2) circuit afflicted by local noise of strength 6^^^; because of the self- 
similarity of our recursive circuit construction, the map from e^^^ to s^^^ is the same as from e to 
e^^\ The noise in the level-(A:— 2) circuit can in turn be coarse grained, thereby level reducing the 
initial level-A: circuit to a level-(A:— 3) circuit. And so on, where at the /-th coarse-graining step the 
input noise strength e^^~^^ is renormalized to an output strength 
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Fig. 1.8. A recursive construction of a multiply encoded quantum circuit. At the base, the unencoded level-0 
circuit corresponds to our quantum algorithm. One level higher, the encoded level- 1 circuit is obtained by 
simulating every elementary operation in the level-0 circuit by using a gadget. By repeating this replacement 
rule, we eventually obtain an encoded level-/c circuit which is the circuit we physically implement. To under- 
stand the effect of local noise on the level- /c circuit, we can level reduce it to a level- (/c—1) circuit which is 
acted by noise of renormalized strength. If level reduction is repeated k times, the level-A: circuit can be level 
reduced to an unencoded level-0 circuit; we can view the strength of noise acting on this level-0 circuit as the 
effective noise strength acting on the multiply encoded operations in the level- /c circuit. 



After k coarse-graining steps, the initial level-fc circuit is eventually level reduced to an unen- 
coded level-0 circuit afflicted by local noise of strength ^(^^ — this level-0 circuit corresponds to the 
quantum algorithm which is simulated by the initial encoded level-A: circuit since every level reduc- 



tion takes us down one level in the ladder in fig. ( |1.8| ). If we use the recursion eqs. ( |1.86| ) where 
^(0) = 5 is the physical noise strength, we find 
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(1.87) 



The constant ^o is the critical noise strength below which the recursive encoding scheme we have 
described is effective; if ^ < sq then s^^^ decreases double exponentially with the coding level k. The 
critical noise strength ^o is often referred to as the threshold for fault-tolerant quantum computation. 
The k successive level reduction steps tell us that we can view the encoded operations in the 
level-fc circuit as being afflicted by an effective local noise of strength £^^\ Thus, the accuracy 1 — 6 



of the level-A: circuit can be estimated as in S ection [L4| where we estimated the accuracy of a level- 1 
circuit afflicted by local noise of strength £^^\ For local Markovian noise, we now have 



6<{e-l)L£^^^ 



(1.88) 



where we used eq. ( 1.85 ) with ^^^^ replaced by £^^\ If the physical noise has a strength below 
the threshold, e <£o, then 6 can become as small as desired by recursively re-encoding the level-0 
circuit of our quantum algorithm sufficiently many times k. The same conclusion also holds for 
local non-Markovian noise. 

In particular, imagine that we desire to obtain the computation output with an accuracy 1 — 6 > 
1 — 6o for some constant (error) Sq, independent of the size L of the quantum algorithm. We can 
arrange to have 6 < Soby choosing k so that 



(t+1)^ < 



log{{e-l)Leo/6o) 
log (eo/e) 



(1.89) 
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Because of the recursiveness of our encoding construction, each of the L encoded operations in 
the level- A: circuit can be implemented by using at most (Lo)^ elementary operations, where Lq 
is the number of elementary operations in the largest gadget. The ratio then of the number L* of 
elementary operations in the entire level- A; circuit over the number L of elementary operations in 
the quantum algorithm scales as 

^ < {L,f = O ilogLf , with a = -i^l^ . (1.90) 

L log(t+l) 

Thus, not only does the level-A: circuit achieve the desired accuracy 1 — 5o. but it does so very 
efficiently; the level-fc circuit is only larger than the unencoded level-0 circuit by a polynomial in 
the logarithm of the size L of the quantum algorithm. 



1.6 Assessment 

The idea in the previous section was to establish the existence of a critical noise strength by con- 
sidering a specific fault tolerance scheme; in particular, we chose to study the recursive scheme 



illustrated in fig. |1.8| because its self-similar nature greatly simplified our analysis. Although recur- 
sive schemes are easier to analyze, it is clearly possible that they are not optimal from a practical 
point of view, and other more complex schemes may have higher thresholds and/or more favorable 
overhead costs. Proposing and analyzing improved schemes for fault-tolerant quantum computation 
is a major focus of current research. 

The existence of a critical noise strength is significant because it implies that the quest to build 
a reliable quantum computer is not a mere fantasy, but it is based on firm foundations: We have 
learned that if we find a physical setting allowing us to experimentally implement quantum circuits 
with local noise of strength ^o or less, then the noisy operations can be assembled efficiently to 
perform an encoded quantum computation and obtain the computation result to as high an accuracy 
as desired. 

Although this knowledge gives us confidence and encouragement to research further how quan- 
tum computers can be constructed, it is possible that the outcome of this endeavor may ultimately 
be failure. We can contemplate several possibilities under which such a failure might occur: First, 
it is possible that the entire concept of what it means to quantum compute, a concept which is based 
on the laws of conventional quantum mechanics, is flawed when applied to quantum computers with 
either a very large number of qubits or very long running times — clearly, in all our considerations 
we have assumed (as the majority of physicists currently believe) that the framework of quantum 
mechanics can be extrapolated without change to the (long) time and (large) length scales rele- 
vant for quantum computers implementing useful computations. Perhaps there are fundamental, as 
yet unknown, principles that prevent the realization of the highly entangled multi-particle quantum 
states required to implement useful quantum algorithms. In this sense, the project of quantum com- 
puting can be seen under a different light, the light of testing quantum mechanics in new regions 
of the parameter space; even if nothing useful as regards computation comes out, we may uncover 
puzzles forcing us to revise our approach to quantum mechanics and physics in general. 

A second possibility for failure relates to the conditions we imposed on the noise as we formu- 
lated our theoretical analysis. It is possible that, as we design and test quantum computing devices 
of increasing complexity, we will eventually find that the physical noise is not captured by our local 
noise models or that, even if noise is local, its strength cannot be upper bounded by a small constant 
number. Ultimately, the question of whether methods of quantum fault tolerance can in practice be 
as effective as our theoretical analysis indicates will be decided by the progress of the future exper- 
iments. In the mean time, theoretical research has still ample room for further progress: Fruitful 
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new research can attempt to relax the requirements under which reliable quantum computation can 
be provably shown to be possible; e.g., one may consider more precise models for the noise during 
qubit preparation and measurement, one may specialize to noise models that more closely describe 
the particular characteristics of observed decoherence in modern prototype experimental devices, 
etc. 

At present, we have no evidence neither that quantum mechanics is violated at the length and 
time scales relevant for long useful quantum computations, nor that the physical noise in prospec- 
tive implementations of quantum computation has features that prevent quantum fault tolerance 
from working. Certainly, there are formidable technical difficulties for building a large-scale quan- 
tum computer with present technology, and it is possible that the engineering requirements may 
prove too challenging to overcome for a long time in the future. Nevertheless, experimental efforts 
during the last decade have shown great progress, and there is a great sense of optimism among 
experimentalists that this progress will continue even more rapidly as they gain more insight and 
intuition about their systems. 



1.7 History and further reading 

Computer engineering is the art and science of translating user requirements we do not fully understand; 
into hardware and software we cannot precisely analyze; to operate in environments we cannot accurately 
predict; all in such a way that the society at large is given no reason to suspect the extent of our ignorance. 
— adapted from Kaplan's By Design: Why There Are No Locks on the Bathroom Doors in the Hotel Louis 
XIV and Other Object Lessons, Fairchild Books (2004). 

In the hope of making the flow of thought in this chapter as smooth as possible, we have avoided 
interruptions to discuss the history of the subject of quantum fault tolerance and we have also omit- 
ted discussing a number of technical but important details. This final section provides some of this 
historical context and references to published works — most of which are available freely on the 
arXiv.org servers — where further information can be obtained. Of course, knowing that our histori- 
cal account and our list of references cannot be perfectly complete, our aspiration is not to provide 
an exhaustive list of all relevant publications but rather to guide the interested reader in his/her first 
steps in the large bibliography. 

The question whether logical operations can be implemented fault tolerantly despite noise was 
central from the early days of the development of classical computing. Shannon's master's thesis [Q 
laid the foundations of digital circuit design, and von Neumann's analysis of noisy cellular automata 
loo showed how unreliable components can be assembled to implement reliable computations. A 
more recent exposition of methods for reliable classical computation can be found in Gacs' work 
PI ; see also Gacs' works [5, 6| on noisy cellular automata and Gray's guide fT] on fS]. 

The corresponding study for quantum computing was pioneered by Shor [ 8 1 who described the 
first gadget constructions for universal quantum computation. Soon after, the existence of a critical 



noise strength based on a recursive scheme as in fig. |1.8| was discussed by Aharonov and Ben-Or 
SUni, by Kitaev lElIEl, and by Knill, Laflamme, and Zurek |[T3l. All these works considered 
local Markovian noise and made a series of additional assumptions about the experimental quantum 
computing devices. Most notably, one assumes that one can supply fresh ancillary qubits or refresh 
existing qubits at any point in time during the noisy quantum computation (this is a necessary 
assumption; if it is dropped, there is no critical noise strength lO). In addition, one assumes 
that there is maximum parallelism, i.e., it is possible to apply gates in parallel on disjoint sets of 
qubits (also a necessary assumption flSl). Finally, one assumes that multi-qubit gates can be applied 
between any set of qubits irrespective of their geometric distance (this is not a necessary assumption; 
a critical noise strength exists even when geometric constraints are taken into account |[T6l ). 
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More recently, a new proof for the existence of a critical noise strength was described by Al- 
iferis, Gottesman, and Preskill [17J; see also Aliferis' doctoral thesis [18] and Gottesman's review 
lfT9ll . This proof is significantly simpler than earlier proofs, and it applies to both Markovian and 
non-Markovian local noise (the analysis for non-Markovian local noise extends prior results by 
Terhal and Burkard EOl ). Building on this new proof, Aharonov, Kitaev, and Preskill ETTl later an- 
alyzed long-range static noise, Aliferis and Terhal Il22l analyzed leakage noise, Aliferis and Preskill 
||23]| analyzed biased noise with dephasing being much more dominant than relaxation or leakage, 
and Ng and Preskill |[24l analyzed Gaussian noise. 

The value of the critical noise strength has been estimated both analytically by means of com- 
binatorial analyses and also, for simple probabilistic local noise models, by performing numerical 
simulations in a classical computer. The highest numerical estimates to date (of order 1.0 x 10~^) 
have been obtained for Knill's postselection and Fibonacci schemes |25 1 and for Raussendorf, Har- 
rington, and Goyal's scheme based on surface codes |[26ll . The highest analytical estimates to date 
(of order 1.0 x 10~^) have been obtained by Reichardt ll27]| , by Aliferis, Gottesman, and Preskill 
ll28ll . and by Aliferis and Preskill 11291 , all by analyzing Knill's schemes and modifications of them. 
It is interesting to note that the schemes with the highest known critical noise strengths share two 
features: First, they make use of quantum teleportation |30| for implementing quantum error correc- 
tion |[25]| and for simulating certain gates lISTIl . Secondly, they use a method by Bravyi and Kitaev 
ll32ll for distilling high accuracy copies of certain ancillary quantum states out of noisier copies of 
the same states. 

- Panos Aliferis {2m9) 
panos @ alumni, caltech. edu 
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